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We analyze the basic hard exclusive processes: 7T7*7-transition, pion and nucleon electromag- 
netic form factors, and discuss the analytic continuation of QCD formulas from the spacelike 
q < to the timelike region q > of the relevant momentum transfers. We describe the con- 
struction of the timelike version of the coupling constant a s . We show that due to the analytic 
continuation of the collinear logarithms each eigenfunction of the evolution equation acquires a 
phase factor and investigate the resulting interference effects which are shown to be very small. 
We found no sources for the if-factor-type enhancements in the perturbative QCD contribution 
to the hadronic form factors. To study the soft part of the pion electromagnetic form factor, 
we use a QCD sum rule inspired model and show that there are non-canceling Sudakov double 
logarithms which result in a /("-factor-type enhancement in the timelike region. 



I. INTRODUCTION 



Within the factorization framework, perturbative QCD has been applied to various processes involving large 
momentum transfers, both in the spacelike q 2 — —Q 2 < (for reviews, we refer to and the timelike 

q 2 > regions (see, for example, [|7H10[). Note that the running coupling constant a s (fjr) is usually defined 
with reference to some Euclidean (spacelike) configuration of momenta of scale /x. For large spacelike q, this 
produces no special complications. One simply uses the renormalization group to resum the logarithmic corrections 
(a s (/j, 2 )\ri(Q 2 / /j, 2 )) N that appear in higher orders of perturbation theory, arriving at an expansion in the effective 
coupling constant a s (Q 2 ) which, in the 1-loop approximation, is given by jjj 

" s(Q ) = (ll-2JV / /3)ln(Q 2 /A 2 ) ' (1) 

with Nf being the number of active flavors and A denoting A^ . In general, the A-parameterization of a s (Q 2 ) 
is a series expansion in 1/L (where L = ln(Q 2 /A 2 )), and the definition of A is fixed only if the 0(l/L 2 )-term is 
added to Eq.(|]) [jllj. Continuing the logarithms into the region of timelike q, one should deal with the in terms: 
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ln(Q 2 //i 2 ) — > ln(Q 2 //i 2 ) ± i7r, which may produce large higher-order corrections. In the case of the i?-ratio for 
e + e~ — > hadrons process, this problem was discussed in refs. []l2[-p^|. It was shown there that, by using the A- 
paramcterization for a s (Q 2 ) in the spacelike region, it is possible to construct for R(q 2 ) an expansion in the timclikc 
region in which all the (tt 2 / L 2 ) N -terms are resummed explicitly, and, what is most important, the transformation 
into the timelike region reduces the magnitude of each particular term of the 1/L expansion. Another well-studied 
example related to the analytic continuation into the timelike region is the cross section of the Drell-Yan (DY) 
process AB — > 7*A. In this case, the iir factors associated with the continuation of the Sudakov double logarithms 
(a s \n 2 (Q 2 /m 2 )) N result in a 7r 2 -enhanced correction which gives rise to the if-factor [[l5| increasing in turn the 
result of the perturbative QCD calculation by the factor of 3 to bring the DY cross section in agreement with 
experiment. 

For elastic form factors, existing experimental data -Q show a considerable enhancement of the timelike form 
factors over their spacelike counterparts. In the present paper, we study possible sources of such an enhancement. 
To disentangle different aspects of the analytic continuation into the timelike region, we proceed step by step, 
beginning with the simplest cases and then going on to more complicated ones. We start with a discussion of 
the analytic continuation into the timelike region of the UV logarithms ln(Q 2 / u^) inducing the Q 2 -dependence 
of the running coupling constant a s (Q 2 ). We take the cleanest case of R{e + e~ — > hadrons), in which no other 
types of logarithms appear and review in Section II the continuation procedure for R(e + e~ — » hadrons) as given 
in refs. (l^Jl^]. In section III, we consider another fundamental process: 7*7 — > ir°. At the leading logarithm 
level, only the collinear logarithms ln(Q 2 / u'p) are important while a s can be treated as a constant. So, this is 
another "clean situation" which gives an opportunity to concentrate on the study of the analytic continuation of 
the collinear logarithms which induce the <5 2 -dependence of the pion distribution amplitude ip n (x, Q 2 ). In Section 
IV, we briefly discuss the effects due to the analytic continuation of the Sudakov double logarithms. We consider 
first the cross section of the Drell-Yan process AB — ► [i + u~X. In this case, the double logs ln 2 (Q 2 //z 2 ) appear on 
a diagram by diagram basis but cancel after resumming over all diagrams of a given order. However, the ir 2 terms 
generated by the analytic continuation survive and, as already mentioned, produce an enhancement due to the 
if-factor. We contrast this outcome with the case of the hard contribution to the pion electromagnetic form factor, 
in which the induced ir 2 terms cancel together with the double logs. For this reason, the modification of the hard 
term of the pion form factor in the timelike region is only affected by the analytic continuation of the UV and 
collinear logarithms. These effects are discussed in Section V. In Section VI, we study the analytic continuation 
of the hard pQCD contribution to the nucleon form factor. Both in the pion and the nucleon case, we find that 
the effects due to the continuation into the timelike region are very small. Experimentally, however, the timelike 
nucleon form factor is essentially larger than its spacelike counterpart. This discrepancy may be regarded as an 
indication that the hard contribution does not dominate the form factors at accessible momentum transfers. An 
alternative scenario discussed in many papers |l9|-^6| is that in the few GeV 2 region the form factors are dominated 
by the soft mechanism. In Section VII, we study the analytic continuation effects for the soft contribution to the 
pion electromagnetic form factor within the local quark-hadron duality model motivated by the QCD sum rule 
analysis of refs. [p^|- ^l|j2^ 1 . We show that at the one loop level, there are explicit non-canceled double logarithms 
\n 2 (Q 2 /fj> 2 ) which produce the n 2 terms in the timelike region, giving rise to a i^-factor-type enhancement. 



II. CONTINUATION OF a s INTO THE TIMELIKE REGION: R(e+e -> hadrons, s) 

The ratio R(s) = a(e + e~ — ► hadrons)/<j(e + e~ — * fi + u~), characterizing the total cross section of e + e~ annihila- 
tion into hadrons, provides the simplest example of the analytic continuation of the effective QCD coupling constant 
a s into the timelike region. The standard procedure (see, e.g., J2^j] and references cited therein) is to calculate the 
Adler function D(Q 2 ) by taking the derivative D(Q 2 ) = Q 2 dH/dQ 2 of the vacuum polarization II(Q 2 ) related to 
R(s) by 

R(s) = ^- (n(- s + ie) - n(- s - ie)) . (2) 

In perturbative QCD, D(Q 2 ) is given by the a s (Q )-expansion: 

In the MS scheme, the coefficients di are known up to i = 3 ]2"7|j28| ]. Using Eq. and the definition of D, one 
can relate i? QCD (s), the perturbative QCD version of R(s), directly to D QCD (Q 2 ) 
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The integration contour in Eq.(^J) goes below the real axis from — s — ie to some point Q 2 in the deep spacelike 
region and then above the real axis to — s + ie, i.e. in the region where the function D(s) is analytic. 

In a shorthand notation, D — > R = $[£)]. The actual calculation is very simple if one represents a s (Q 2 ) through 
an expansion in 1/ ln(Q 2 /A 2 ), i.e., via the A-parameterization. The latter results from the QCD Gell-Mann-Low 
equation 



L = ln(Q 2 /A 2 ) 
where bk are /^-function coefficients: 
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6 = 11 — 2JV//3 @|, 6i = 102 - 38iV//3 ||, bf s = 2857/2 - 5033^/18 + 325iV 2 /54 @. 

Inverting (|5|) by iterations and reexpanding the result in 1/L we get the A-parameterization for the running 
coupling constant 
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where L\ = (6i/6q) l n (^o-^) — A pl| . ^2] |. To fix the functional dependence of a s {Q 2 ) on Q 2 , one should specify the 
integration constant A. The standard (or "popular") choice is 



A po P = ln5o 
K 



(7) 



which gives the shortest expression ln(L) for L\, A clear disadvantage of this choice is that it guarantees a 

rather large 1/L 2 correction to a s , which results in a large difference between A LO and A NL °. As argued in Rcf. 
[fl3[ , a more appropriate (optimal) choice is 



A°p* = ^ln6 Z 



(8) 



where L is the average value of the logarithm L within the region under study, e.g., L = 4 corresponding to 
a s /ir ~ 0.1. For this choice, the ratio Li/L is smaller than 7% and Eq.( ||) has 1% accuracy in the whole region 
L > 3, with the total correction to the simplest formula (1) being less than 10%. 
The A-parameters corresponding to different A's are related by 
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In particular, 



A opt = A pop /L bl / 2b o 



(9) 



(10) 



Taking L = 4 we get A op '|^ =4 w A pop /1.73. In connection with the discussion above, we want to stress here that 
preparing to analytically continue an approximate expression it makes sense to take care of the convergence quality 
of the original expansion in the spacelikc region. If there are corrections which are under our full control and we 
can make them small, then we should use this opportunity and make them small. 

Now one can substitute a s (Q 2 ) in Eq.(||) by its A-parameterization to get an 1/L expansion for the Adler function 
D(Q 2 ). For each term of this expansion, the integral ([IJ) can be calculated explicitly (see also |53j) 



1-1, 

► — (n/2 — arctan(L s /7r)) 

L„ TT 



i>A 2 



If 1 TX Z 

- arctan(vr/ J L s ) = — <j 1 - - — 



(11) 
(12) 



ln(L CT /L ) ln(V i 2 . + tt 2 /L ) - (L s /ir){n/2 - arctan(L s /7r)) + 1 



LI 



L 2 



i>A 2 



ln(^/L 2 + 7r 2 /L ) - (L s /n) arctan(7r/i ;j ) + 1 _ L s /L 



L 2 



LI 



7T 

I! 



O TT 

6IA 



+ 



(13) 



3 



1 1 1 



1-TO+--- , (14) 



jji ip, _|_ n 2 L 2 \ L 2 

Z» ~* ^ (n-1)! (dZj) Zl + 7T 2 = I" j 1 ~ If 6 + "'} ' ^ 

where L s = ln(s/A 2 ), L a = ln(er/A 2 ), and we assume that s > 0. Furthermore, Lq = e Ab °/ bl /bo is a constant 
depending on the A-choice in the A-parameterization. 

Using (g) and incorporating Eqs.([il"|)-(|l5|) (as well as their generalizations for In 2 L/L 3 , InL/L 3 , etc.) one obtains 
the expansion for i?^ CD (s) 

R QCB {s) = £ e 2 j 1 + £ ^[(^/^ I (16) 
g I fe=i J 

in which all the (7r 2 /L 2 ) jv -terms are resummed. 

As noted in Ref. |13[, the application of the ^-operation normally violates nonlinear relations: $[1/L 2 ] ^= 
($[1/L]) 2 , etc. However, it respects linear relations $[A + B] = $[A] + $[£?], $[\A] = A$L4] and 
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In particular, this relation was used to explicitly obtain $[1/X n ] in Eq.([l5|). As a result, expansion ( |l6| ) is not 
an expansion in powers of some particular parameter. A priori, there is no reason to believe that a power series 
expansion is better than any other. In fact, expansion (|l^) converges better than the generating expansion (Q) for 
D(a) because, as it follows from Eqs. (p^-(|l5|), $[af] is always smaller than . Moreover, ($[a^ +1 ]) 1 ^ Ar+1 ^ < 
($[af })^ N , i.e., the effective expansion parameter decreases in higher orders. Thus, if one succeeded in obtaining 
a good expansion for D(a) (with all cIn being small numbers), then the resulting $[a^] expansion for i?Q CD (s) 
is even better, and the best thing to do is to leave it as it is. 

The timelike analogue of the simplest A-parameterization for a s (Q 2 ) (Eq.(Q)) is then 
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This function has a finite value both at q 2 = A 2 and q 2 = 0. 

The well-known deficiency of the perturbative expansion for DQ CD (Q 2 ) in powers of a s (Q 2 ) is the presence of 
the unphysical singularity at Q 2 = A 2 induced by the Landau pole of 1/ ln(Q 2 /A 2 ). As a consequence, 
as calculated from Eq. (^), also has unphysical features: namely, it does not vanish on the negative real axis. In 
particular, substituting 1/L a into the integral (§) and taking negative s we get 



1 



(-A 2 < s < 0) , (19) 
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which results in an unphysical cut of II ( ^ CD (s) in the region — A 2 < s < 0. Furthermore, applying Eq. (^) to the 
pole term £P ole (Q 2 ) = A 2 /(Q 2 - A 2 ) one obtains the result coinciding with the rhs of Eq.([l9|). Hence, if one now 
postulates that D^ CD (<5 2 ) is given by integrating i?Q CD (s) over the physical region s > only, i.e., if one takes 



(20) 



(this transformation will be denoted as R — > D), then £)Q CD ((3 2 ) is free from the unphysical singularities at 
Q 2 = A 2 . For instance, combining the two transformations (D — ► R — > D) = (D =>• D) one would get 
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which coincides with the pole-free expression for the running coupling constant proposed by Shirkov and Solovtsov 
However, since the D — > R operation does not respect nonlinear relations, the D D transformation acting 
on would not produce the nth power of the rhs of Eq. (pi|) . Hence, a s cannot serve as an expansion parameter 
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of a power series. Noting that both D — > R and R — > I) convert derivatives with respect to the logarithm of the 
initial variable into derivatives with respect to the logarithm of the resulting variable we obtain 

1 1 d n ~ 1 1 1 d n ~ x { 1 A 2 \ 

_ = (-1)" ^ r— IV" - - I— - I (22) 
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This relation was given in a recent paper by Shirkov f35| , see also Ref. p6| for a related discussion of perturbation 
theory expansions in the timelike and spacelike regions. 

For moderate values of Q 2 , the modification due to the continuation into the timelike region is numerically 
rather significant: for a s > 0.3 the 7r 2 /L 2 -terms change a s by more than 20%, i.e., they are more important (for an 
optimal choice of the A-parameter) than the 2-loop corrections in the A-parameterization (|^) . On the other hand, 
the difference between d s (Q 2 ) and the modified spacelike coupling u s (Q 2 ) (taken at mirror momenta) is rather 
small (less than 10%) for all values of Q 2 . 

Thus, using the A-parameterization for the effective QCD coupling constant in the spacelike region, we obtained 
an explicit expansion for the timelike quantity RQ^(s). One may question, though, the reliability of the above 
formulas in the region of small momenta \q\ ~ A. In particular, a rapid change of a s in the small-g 2 region (compare 
a s (A 2 ) — 27r/6o and a s (0) = Air/bo) is as suspicious as the Landau pole of a s (Q 2 ). Evidently, they both are artifacts 
of the analytic continuation procedure applied outside the applicability region. It is well known that the physical 
R(s) vanishes below the two-pion threshold and approaches the perturbative value only for values of s marginally 
larger than A 2 . So, one may argue that a more realistic procedure is to integrate R^ CD (s) in the dispersion relation 
( |20|) from some effective threshold sq rather than from zero. Taking, e.g., sq = A 2 , one would get another effective 
spacelike coupling, call it a s (Q 2 ). It vanishes at Q 2 = 0, but is essentially constant a s {Q 2 )/ii 0.1 in a wide range 



A 2 < Q 2 < 30A 2 of momenta. Hence, a s (Q 2 ) effectively "freezes" at small momenta (see also [B7] B9 33 ). 



III. COLLINEAR LOGARITHMS AND DISTRIBUTION AMPLITUDES IN THE TIMELIKE REGION 

The logarithmic dependence on the large momentum scale Q 2 may also appear through mass logarithms 
ln(<5 2 /m 2 ), where m is some mass or an infrared regularization parameter. Note that the standard pQCD factor- 
ization 

T(Q 2 /m 2 ) = t(Q 2 /fi 2 ) (g) ipifj 2 ) (23) 

works only in a single-logarithm situation, when there may appear just one ln(Q 2 /m 2 ) factor per each loop. These 
collinear logarithms can be absorbed into the renormalization of the long-distance function (distribution amplitude) 
tp([i 2 ). In particular, taking /i 2 = Q 2 , one arrives at the description in terms of Q 2 -dependent functions (p(Q 2 ). 
Again, if the large momentum is timclike, the collinear logarithms ln(Q 2 /m 2 ) acquire the imaginary part ±i7r, and 
we may ask how one should define the Q 2 -dependent distribution amplitudes f{Q 2 ) in the timelikc region. 

To approach this problem, let us consider the simplest example of a hard exclusive process: ir° production in 
7*7 collisions. Its pQCD expansion starts at zero order in a s 

t (x,Q 2 ) = -y, (24) 
and the leading pQCD result [[l0| for the large-Q 2 behavior of the form factor is 

TP frflS 47r f 1 <P*i X ) J - 47r A r / rN 

^ (Q) = li ~x~Q r = ~3Q 2 (25) 
The nonperturbative information here is accumulated in the same integral 
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that appears in the one-gluon-exchange diagram for the pion electromagnetic form factor |41|~|43|. The value of / 
depends on the shape of the pion distribution amplitude tp^ix). In particular, using the asymptotic form |p~l| , fi"2"f 

yf{x)=&Ux{l-x) (27) 

gives /q S = 3. If one takes instead the Chernyak-Zhitnitsky model Q] 
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<pf{x) =30/^(1 -x)(l- 2a;) 2 , 



(28) 



the integral Jo increases by a sizable factor of 5/3: Iq Z = 5. This difference can be used for an experimental 
discrimination between the two competing models for the pion distribution amplitude. 

At one loop, the MS coefficient function for the 7*7 — > n° form factor was calculated in refs. fl45|-[47t and was 
found to be 
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In full compliance with the factorization theorems [|48 40 (see also [[I9|-pl|), the one- loop contribution contains 
no Sudakov double logarithms In 2 Q 2 of the large momentum transfer Q. Physically, this result is due to the color 
neutrality of the pion. In the axial gauge, the Sudakov double logarithms appear in the box diagram but they are 
canceled by similar terms from the quark self-energy corrections. In Feynman gauge, the double logarithms In 2 Q 2 
do not appear in any one-loop diagram. It is easy to check that the term containing the logarithm ln(Q 2 / ^ 2 ) has 
the form of a convolution 
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of the lowest-order ("Born") term to(^,Q 2 ) = 1/£,Q 2 and the kernel 
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governing the evolution of the pion distribution amplitude. The "+" -operation is defined here, as usual 



(31) 
by 



[F& x)]+ = x) - 8{i - x) J F(C, x) dC . 
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Since the asymptotic distribution amplitude is the eigenfunction of the evolution kernel V(£, x) corresponding 
to zero eigenvalue 



V(£,x) <p as {x)dx = 0, 



(33) 



the coefficient (| + lnx) of the ln(Q 2 //i 2 ) term vanishes after the ^-integration with <^ as (x). Hence, the size of 
the one-loop correction for the asymptotic distribution amplitude is /^-independent and is determined only by the 
remaining terms (for a detailed discussion of their structure, see Ref. 

In this section, we want to concentrate on the Q 2 -dependence induced by collinear logarithms, which in this 
process start to appear at the one-loop level. The UV logarithms shifting the argument of a s appear only at 
two-loop order. Hence, analyzing the leading collinear logarithms (a s ln(Q 2 / fi 2 )) N we will treat a s as a constant. 
The factorization theorem means essentially that the leading logarithms (a s \n(Q 2 / fi 2 )) N exponentiate in higher 
orders producing a factor which can be absorbed into the renormalization of the pion distribution amplitude 
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Now, taking a timelike momentum Q 2 = —q 2 , we would get an extra ±iir term: ln(Q 2 //^ 2 ) — * ln(q 2 //j 2 ) ± iir and 



exp[- ln(Q 2 /V 2 ) V] -> exp[- ln(g 2 /^ 2 ) V] exp[±iirV) . 



(35) 



The first exponential corresponds to the standard evolution of the pion distribution amplitude from the scale [i 2 
to the scale q 2 . The second exponential is specific for the timelike kinematics. In our approximation, it is q 2 - 
independent and can be treated as a conversion factor for the transition from a "spacelike" distribution amplitude 
(p to its timelike counterparts (p± 



<p± = exp[±«7ry] (8 f ■ 

In general, "timelikc" distribution amplitudes have both real and imaginary parts. However, since V ® ip' 
the spacelike asymptotic distribution amplitude does not differ from its timelike counterpart. 
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To estimate the effect of phases, let us consider the case when the spacelike distribution amplitude is given by 
the Chernyak-Zhitnitsky (CZ) model 0], which can be represented as 
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where <^ as = 6/ ff x(l — x) and P2 = 24/ 7r a;(l — x)(l — 5x(l — x)) is the next eigenfunction of the V kernel 
corresponding to the eigenvalue 72 = ?§a: s /7r. The timelike distribution amplitude is then 



^CZ = ^as + e ±i(25/18)a. ¥ , 2 
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and the / integral for this function is 



l9 z = 3 + 2e ±l ^/^ c 
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Its absolute magnitude 
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is slightly smaller (by 2% if a s — 0.3) than the spacelike value I 



IV. SUDAKOV LOGARITHMS AND if -FACTOR 



Small radiative corrections in the timelike version of the 7*7 — > 7r° process are in strong contrast with the large 
-ftT-factor value found for the Drell-Yan process AB — > 7*Y. These corrections originate from the Sudakov double 
logarithms (a s ln 2 (Q 2 /m 2 )) N . In the spacelike region, the double logarithms due to the virtual gluon exchanges 
exponentiate into the Sudakov form factor 

S(Q 2 /m 2 ) = e- a ° i" 2 W 2 /™ 2 )/3- (41) 

(again, we treat a s as a constant). In the DY process, the photon momentum is timelike, and the logarithm 
ln(Q 2 /m 2 ) acquires the ±iir additional term, so that one has 

- L 2 -> -L 2 ± 2iirL + 7T 2 . (42) 

The imaginary parts of the two conjugate diagrams shown in Fig.|l|a,b cancel, the double log L 2 from Fig.|l|a 
(b) is also canceled by the real gluon emission diagram Figjljc (d), while the 7r 2 -term survives and leads, after 
exponentiation, to a large K factor exp[27ra s /3] ~ 2. The crucial technical observation here is that the real 
emission diagrams give L 2 without 7r 2 -terms. This can be easily understood looking at the reduced diagrams for 
the virtual vertex correction and real gluon emission. Take for definiteness, the Feynman gauge. Then the virtual 
vertex correction diagram Fig.|l|a contains the — ln 2 (— s/m 2 ) term, where s — (xpa +UPb) 2 — xyS is timelike, and 
the resulting contribution contains a tt 2 term. The real emission diagram Fig.[lp, in turn, contains the ln 2 (— u/m 2 ) 
term, where u = (xpa — VPb) 2 = ~xyS is now spacelike, and there is no it 2 term in this contribution. 




FIG. 1. One- loop reduced diagrams for the DY process cross section calculated as imaginary part of the forward scattering 
amplitude, o, b) Virtual vertex corrections, c, d) Real gluon emission corrections. 

For the hard pQCD contribution to the pion electromagnetic form factor (which is considered in more detail 
in the next section), the situation is completely different. In this case, the initial state of the hard subprocess is 
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represented by a qq pair with momenta xp and (1 — x)p. After the hard scattering subprocess, one deals with 
a qq pair with the final momentum p' shared in fractions yp' and (1 — y)p' '. In Feynman gauge, the double 
logarithms \n 2 (Q 2 //i 2 ), where Q 2 = — (p — p') 2 , appear when the reduced diagrams have the structure of those 
shown in Fig.|2j. One can easily check that the relevant momentum transfers in all four cases have the structure 
tij = {xip — VjP 1 ) 2 = —XiUjQ 2 , resulting in the double logs In (— iy//z 2 ). When the momentum transfer q = p' — p 
is spacelike, all iy's are spacelike, whereas for a timelike q, all ty's are timelike as well. In the latter case, one 



has 7T terms for each particular diagram. The double logarithms in diagrams 
because the soft gluon interacts in the final state with quarks of opposite co 

neutrality of the pion, the double logs ln 2 (Q 2 //i 2 ) cancel for the sum of the diagrams of a given order. For timelike 
q, they cancel together with the accompanying it 2 terms. 



Ja and ||b (0c and ||d) differ in sign 
or charge. Hence, due to the color 




FIG. 2. One-loop reduced diagrams for the hard pQCD contribution to the pion electromagnetic form factor, which 
contain the Sudakov double logarithms in Feynman gauge. 

Thus, even for a timelike momentum transfer q, there is no if-factor for the pQCD hard contribution to the pion 
electromagnetic form factor. After cancellation of the Sudakov double logs, only the evolution-related collincar 
logarithms remain, and the situation is rather similar to the simplest case of the 7*7 — * 7r° form factor. 



V. PION FORM FACTOR IN THE PERTURB ATIVE QCD APPROACH 



The general pQCD factorization formula for the pion electromagnetic form factor at large momentum transfer 
reads 
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dx / ip{x, n F ,n R ) T(x, y, Q , fj, F , p, R ) ip(y, fi F , fi R ) dy , 



(43) 



where \i F is the factorization scale for the collinear logarithms and (i R is the renormalization scale for the UV 
logarithms. The hard scattering amplitude is given by an expansion in a s 
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2TTC F a s (p 2 R ) 
xyQ 2 N c 



1 + ^M) T ( 1)(:C!y)Q 2 )M | )M 2 )+o(a 2 



(44) 



The one-loop correction TW(i, y 1 Q 2 , p, 2 F , p? R ) was calculated using the dimensional regularization in several 
papers [|54|-[59| which differ from each other by a particular choice of renormalization and factorization prescriptions. 
These differences (and also typos and mistakes) were discussed in Refs. p^- |60| , p9[ | . In the MS subtraction scheme, 
supplemented by the requirement that both the a s and tp-nix) are process- independent functions, the one-loop 
correction has the form 



T« = C F T F (x, y, Q 2 , M |) + I TP( X , y, Q 2 , + (C F - N c /2) T A (x, y) 



T 1 



3 + ln(xy) 



In 



UJ + 2 ln {xy) + 2 Hxy) 2(1^0 " W^V) 
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~3~ 



T0 = -\n[\)-\n{xy) + °- , 



I'r 

T A = Li 2 (l - x) - Li 2 (x) + ln(l - x) ln 
1 



(45) 



1-1/ 



(x - y) 2 



(x + y — 2xy) ln(l — x) + 2xyhi{x) + 



(1 - x)x 2 + (1 - y)y 2 



x-y 



8 



ln(l - x) \n(y) - Li 2 (l - x) + Li 2 (x) 



+ {x^y} (46) 



(we use here notations similar to those of Ref. |61| ), As usual, Li 2 is the dilogarithm (Spence) function. 

We already discussed in the previous section that all the Sudakov double logarithms ln 2 (Q 2 //x 2 ) cancel and that 
only the collincar single-logarithms ln(Q 2 //ip) remain. Comparing Eq. ( [i"6|) with the one loop correction to the 
7*7 — * 7T° hard scattering amplitude, Eq. (|29|), one can easily notice many similarities in the structure of the 
one-loop corrections in these two cases. In particular, the coefficient [3 + ln(xy)] in front of the evolution logarithm 
ln(Q 2 / ' n 2 F ) is the sum [3/2 + In a;] + [3/2 + In y] of terms corresponding to the convolution of the tree level term 
with the kernels V(x, £)5(r/ — y) and S(x — $,)V(r), y) (see Eq. ( |30| ) ) inducing the evolution of the pion distribution 
amplitudes tp(x,n%) and <p(y,fXp), In a sense, the collinear logarithms indicate that the pion structure is probed 
at a scale proportional to Q. However, one should remember that since the asymptotic wave function does not 
evolve, the coefficient accompanying the evolution logarithm ln(Q 2 / fip) vanishes if the pion wave function has the 
asymptotic shape. As a result, the choice of fiF in that case does not affect the size of the one-loop correction. 
The latter comes from several sources which can be identified in a way similar to the detailed analysis of the one 
loop correction for the 7*7 — ► 7r° form factor given in Ref. [53]. 

In addition to the evolution term proportional to \x\{Q^]x 2 F ), there is a rather large positive correction due 
to the ^ln 2 (xy) term and even larger negative contributions corresponding to the constant term —14/3 and the 
logarithmic term I \n(xy). As explained in Ref. M], in the 7*7 — > n° case, the | In 2 x term is a result of a positive 
In 2 x evolution-related contribution and a negative — | In 2 x Sudakov-related term. As we emphasized earlier, the 
Sudakov In 2 Q 2 double logs should cancel, otherwise there is no pQCD factorization. However, when several scales 
are involved, like Q 2 and xyQ 2 in our case, there may be a remnant like ln 2 (xy). In the pion form factor case, there 
is another scale xQ 2 , the quark virtuality, whence the single logarithms In a; +lny. The latter give a rather large 
negative contribution. There are also large negative constants (—9/2 in the 7*7 — * tt° case and —14/3 in the pion 
form factor case), which are another (and numerically very important) manifestation of the Sudakov effects in the 
impact parameter space. In full analogy with the results of Ref. ^3|, these [and the In 2 (xy),\n(xy) terms] result 
from convoluting the &-space version Ko(y/ ' xyQ 2 b 2 ) of the one-gluon exchange propagator and the &-space Sudakov 
form factors S(x, bQ), S(y, bQ) (exact one-loop expressions are given in |53|]). In the practically important case of 
the asymptotic wave function, the total correction due to the T F term is negative and equal to — (71/18)a s /7r; as 
one could expect, it is approximately twice larger than that in the 7*7 — > tt° case. 

The situation is reversed in the case of the UV related T' 3 term: it is dominated by large positive contributions. 
In full accordance with the renormalization group, the UV logarithm h\{Q 2 / ^ 2 R ) is accompanied by the /3-function 
coefficient bo . It generates the running of the effective QCD coupling a s , "suggesting" that we should use some 
scale proportional to Q 2 as its argument. According to Brodsky, Lepage and Mackenzie |6^], one should choose the 
argument of the effective coupling constant in such a way as to absorb all the terms proportional to bo from the 
next loop correction. Taken literally, the BLM prescription in our case corresponds to using the (rescaled) gluon 
virtuality 

/4 = xyQ 2 e" 5 / 3 « xy(Q/2.3) 2 (47) 

as the argument of a s . The rescaling factor e 5 / 6 w 2.3 reflects the fact that the MS scheme measures the momenta 
in "wrong" units. To cure this effect, one may introduce a version of the minimal subtraction scheme which 
measures momenta in more "physical" units ApjjYS — e5|,6 ^MS' This choice is similar to using the ay coupling of 



Brodsky et al. J39J (note that their relation a v (Q) = a* /s (e" 5/6 Q)[l - ^N c af s /it] includes a NLO correction). 
One should remember, however, that the actual expansion parameter for switching from the leading to the next-to- 
leading level is 1/ ln(Q 2 /A 2 ) rather than a s as a whole. As a result, the "non-physical" nature of the M S scheme is 
almost totally compensated by the non-optimal "popular" choice for the analytic form of a s (Q 2 ). As discussed in 
Section II, A opt w A pop /1.74. As a result, Ap^ Y s ~ ^-^^j^- Due to tne compensation of two opposite corrections, 
the standard A^ji parameter is rather close to the genuine "Aqcd" parameter of the PHYS opt scheme in which 

the coupling a s (k 2 ) corresponding to the vertex with the gluon momentum k is given by 4-7r/&o ln(fc 2 /A 2 ) without 
sizable next-to-leading order corrections. In other words, using the NLO expression for a s in the popular form is 
equivalent to adding a negative term — (bi/bfyhiL to T' 3 , partially compensating the "5/3" constant. For L ~ 4, 
this reduces 5/3 by a factor of 3. Choosing "PHYS" vs. MS and "opt" vs. "pop" one reduces both types of 
corrections which iterate in higher orders. As stated earlier, if the size of some corrections is under our control, it 
is preferable to keep them small rather than rely on cancellation of large terms. The closeness of Ap^ Y s *° ^wj} 

means that discussing the pQCD applicability region one should compare the A 2 parameter of the MS V ° P scheme 
with the actual (unrescaled) gluon virtuality xyQ 2 . 
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However, taking the argument of the effective coupling constant proportional to xyQ 2 one faces the following 
problem: since the integration is over all the momentum fractions in the range < x,y < 1, the "short-distance" 
amplitude in this case always gets contributions from the infrared region of arbitrarily small virtualities. In this 
sense, such an "inside the integral" BLM prescription contradicts the spirit of the pQCD factorization ideology 
which aims at a perfect separation of the short-distance and long-distance effects (at least in perturbation theory). 
The consistent pQCD approach is to apply the BLM prescription to the form factor as a whole, i.e., "outside the 
integral". In this case, one should choose [ir from the requirement that one should get zero for 

-ln(Q 2 /ti)-Hxy) + ^ 

where the "averaging" procedure (...) stands for integration with ipT T (x)ip 7r (y)/xy. This gives a universal x,y- 
independent scale fiR — cirQ, which depends now on the shape of the distribution amplitude. For the asymptotic 
wave function, as we have seen, the average value of Inx is —3/2, hence the "outside the integral" BLM scale is 
(see also [|§) 

^|^as = Q 2 e- 5 / 3 - 3 « (e- 5 / 3 Q 2 )/20 . (48) 

As argued above, in the "pop" treatment, the factor e~ 5 / 3 is largely compensated by the NLO corrections to 
a s (Q 2 ), and, hence the essential virtuality of the "hard gluon" exchanged between the quarks is "only" by a factor 
of 20 smaller than Q 2 , the nominal momentum transfer to the pion. Nevertheless, despite this sizable rescaling 
factor, the pQCD factorization approach is fully consistent in the asymptotic sense: for a sufficiently large Q 2 
one can calculate the short-distance amplitude perturbatively in terms of an arbitrarily small expansion parameter 
a s (Q 2 j '20). For comparison, in the case of Chernyak-Zhitnitsky wave function 

7 

(\nx)\ v cz =-- 

and a^ z w a^ s /2.3: the essential gluon virtualities are 100 times smaller than Q 2 . In this case, one should not 
expect early applicability of pQCD. 

We would like to emphasize that the reason for such a drastic shift of the BLM scale to very low /j,r values is the 
positive large value of the T@ correction: for = Q and ip^ = <^ s , the T 13 term contributes the NLO correction 
10a s /7r. One may be tempted to combine the large positive T@ term and a sizable negative T F term to end up with 
a smaller total correction ~ 6a s /7r. Physically, though, these corrections have a completely different nature: as 
argued above, the T F term comes primarily from the Sudakov effects. Since the latter exponentiate, one deals here 
with the e~ kpas type series in which the sign of the corrections alternates. On the other hand, the UV corrections 
form a geometric series summed into 1/(1 — kpa s ). Hence, there is no doubt that a partial cancellation of the a s 
terms will be followed by an amplified total correction at the a 2 level. Leaving the physically unrelated Sudakov 
and UV corrections separate and addressing the region of experimentally accessible values of Q 2 < lOGeV 2 , one 
should take a s at an infrared scale ~ A (where it freezes at a value close to 0.3) and supplement the result by the 
exponential ~ e ~ 4as / 7T ~ 0.7 of the negative one-loop corrections induced by the 6-space Sudakov effects. Turning 
to the timelike momenta, we cannot find any sources of enhancement: for the |A;| ~ A region we see no other choice 
rather than to take the frozen value a s ~ 0.3 both in the spacelike and timelike regions, while for large momenta, 
the continuation of a s {a 2 F Q 2 ) converts l/L into 1/(L + in) and the ratio {F^f^l/I^t^ is l/V 1 + 7r V- L2 > 
i.e., the timelike term is suppressed compared to the spacelike one. 

Since the structure of the evolution corrections for the pion form factors is essentially identical to that of the 
7*7 — > 7T° form factor, to continue the evolution logarithms, we can use the approach outlined in Section III. The 
first step is to write the solution of the evolution equation [§i|fll 40 as an expansion over Gegenbauer polynomials 
Cn^ 2 (2x — 1), the eigenfunctions of the LO kernel (x, y; a s ): 



(x,Q 2 ) = ^(x) 



fc=i 



(49) 



where as (a;) = 6x(l — x) is the asymptotic distribution amplitude of the pion. The Gegenbauer moments a^iQ 2 ) 
have a simple evolution pattern 

4 fc (Q 2 ) = 4fe(Mo) exp [-72* In (Q 2 /Mo)] (50) 

(we treat a s as a constant here), with j2k being the corresponding anomalous dimensions and a2fe(/ig) the Gegen- 
bauer moments of the initial distribution amplitude 
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l + E a 2fc(Mo)cA(2z-l) 



fe=i 



(51) 



This representation is very convenient to perform the analytic continuation to the timelike region of Q 2 . Indeed, 
changing Q 2 — > — q 2 , one obtains the natural shift ln(Q 2 ) — > ln(g 2 ) + in, so that 



where 



al k {-q 2 ) = aZ k (\q 2 \)e-^ 



(52) 



82k = 7T72fe 



(53) 



From (|52j) it is obvious that the only interesting and potentially enhancing effect is due to the phases 82k, since 
they can destroy some fine tuning of the coefficients a^feC^o) anc ^ produce a positive interference. But in order to 
realize this possibility, one should start with a situation when there are negative coefficients, say, a\ (/i 2 ) < 0, while 
the corresponding phase is close to tt, e.g., 82 ~ 7r. Such a situation is hard to imagine. Even the (unrealistic) CZ 
distribution amplitude has a\ (/i 2 , = (0.5 GeV) 2 ) = 2/3, while other models are closer to the asymptotic distribution 
amplitude, though all models provide et2 > 0. Furthermore, the phase 82 is 25/18a s , so one needs a prohibitively 
large value a s ~ 2.5 for the coupling constant. In Fig. [| we plot the ratio TL/SL for the pion form factor in the 
CZ model, taking the frozen value a s (Q 2 ) = 0.3. As one can see, the absolute value of F^ ARD (g 2 ) in the timelike 
region is reduced. 
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The ratio of the timelike form factor \F. 
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-t)\ over the spacelike one 



case of the CZ model (dashed line) for the pion distribution amplitude 
with a s (Q 2 ) fixed to 0.3. The solid line represents the result for the asymptotic distribution am- 
plitude. 



To conclude, the perturbative contribution to the pion form factor, F^^ RD (Q 2 ), with a realistic distribution 
amplitude (which is close to the asymptotic one, see p3|-|66[) produces no sizable effects in analytically continuing 
to timelike q 2 values. The only potential effect is due to the substitution a s (Q 2 ) — > a s (q 2 ) which results in a 
10%-reduction of the form factor. 



VI. EVOLUTION PHASES OF THE NUCLEON DISTRIBUTION AMPLITUDES AND THE 

NUCLEON FORM FACTOR 



The nucleon form factor in the leading a s order can be cast in the form 

7 



Q 4 GS(Q 2 ) - ^ [4na s (Q 2 )] 2 \h\ 2 J\dx] j\ 



i=l 



(54) 
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where the amplitudes Ti (xj , yj 

) = ^ N) {{x},Q 2 )T^{{x},{y})^ N \{y},Q 2 ) represent convolutions of with the 
appropriate distribution amplitudes cf>^ N \{x}, Q 2 ) evaluated term by term for each contributing diagram (marked 
by the index "i"). The nucleon distribution amplitude can be represented as an expansion over symmetrized 
combinations $({a;}) of Appell polynomials (for more details, we refer to Pf37|]) 



with 



and eigenfunctions 



^ N H{x},Q 2 ) = ^({x})J2B n (Q 2 )<i> n ({x}) 



n=0 



B n (Q 2 )=B n (^ )exp [-7„ln (Q 2 /a«o)] 



rn+n— M 



(55) 



(56) 



(57) 



where J- mn (5, 2, 2; x\, X3) are the Appell polynomials [[]. Here 

'3 



%{M) 



4tt V 2 



C F + 2r jn (M)C l 



(58) 



are the associated anomalous dimensions of trilincar quark operators with the quantum numbers of the nucleon 
containing external derivatives |57|]. In Eq. ( |55| ) as ({x}) = 120 X1X2X3 denotes the asymptotic distribution 
amplitude of the nucleon and B n (iJ^) are expansion coefficients for some initial distribution amplitude 



(59) 



n=0 



Again, the representation given by Eq. (|5^) is very convenient to analyze the analytic continuation of the (hard 
part of) the nucleon form factor into the timelike region of Q 2 . Continuing Q 2 — > — q 2 one obtains in Eq. ( |56| ) 
the same shift as in Eq. (52) with 5 n = a s ([3o/ A)"/ n . Specifying the particular values of the coefficients B n we can 
calculate the ratio of timelike to the spacelike form factors for several models known in the literature. In Fig. |i| we 
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FIG. 4. Ratio of the timelike (with t = q 2 ) to the spacelike (with t = Q 2 ) hard form factors of the nucleon using three 
different nucleon distribution amplitudes (Chernyak-Ogloblin-Zhitnitsky, heterotic and Gari-Stefanis models). 



display the ratio of the timelike to the spacelike form factors of the nucleon for three different nucleon distribution 
amplitudes: Chernyak-Ogloblin-Zhitnitsky, heterotic and Gari-Stefanis |6[]69]|. From this figure, it is clear that 
there is no enhancement due to the analytic continuation except from a marginal factor of order 1.03 for the GS 
model. 



*One can also expand &n({x}) over the polynomials proposed in Ref. ra|; the particular choice of the basis is not essential 
for our purposes. 
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VII. SOFT TERMS FOR THE PION FORM FACTOR IN THE LOCAL DUALITY APPROACH 



So far, we have discussed only the hard pQCD contributions to the hadronic form factors. But as argued in Rcfs. 
|L9-25|, the dominant contribution at intermediate values of the momentum transfers Q 2 < 10 GeV 2 is generated 
by the soft contribution which involves no hard gluon exchanges. As a model for the soft contribution, we assume 
the local duality (LD) approximation in which it is assumed that the pion form factor is dual to the free quark 
spectral density [ |2l| JzOj ] 

f 2 F„(Q 2 ) = i / ° f ° p (s,s',Q 2 ) ds ds' . (60) 

T JO Jo 

The latter is given by 



p (s,s',t) = ~ 



dt 3 \dt 



- yj(s + s' + t) 



1 Ass' 



(61) 



Here the duality interval so corresponds to the effective threshold for the higher excited states and the "continuum" 
in the channels with the axial current quantum numbers. 

In principle, the value of so is fixed by the ratio of the nonperturbative power corrections to the (leading) 
perturbative term in the OPE for the correlator. In what follows, we use the value so ~ 0.7GeV 2 which has been 
extracted in the pioneering paper Q] from the QCD sum rule analysis of the correlator of two axial currents. The 
LD prescription for this correlator just implies the relation 

so = 4tt 2 /' , (62) 
between Sq and the pion decay constant / w . This relation ensures that the Ward identity for the pion form factor 

F n (0) = 1 (63) 

is satisfied within the LD approach. 

Performing the integral on the rhs of Eq. ( |60| ) we get the explicit analytic expression for the pion form factor 

F LD/ n 2x _ , l + 6s /Q 2 

originally obtained in pj| , fr0| . Note that for t > 0.6 GeV 2 , expression ( pi}) is in good agreement with existing data 
(see Fig. |). 

A simplified version of the LD model is based on using the "duality triangle" instead of the "duality square" . In 
this approach [|l9|,|o|,|6]] , one uses the variables S = s\ + s 2 and si — S2, introducing the reduced spectral density 

,s 

p (S,Q 2 )= / Po (S-s',s',Q 2 ) ds' . (65) 



The LD relation, eq. (60), is then substituted by its "triangle" version (TrLD) 

F^Q 2 ) ~ ^ TrLD (Q 2 ) = ^- [ S ° po(S, Q 2 ) dS (66) 

^ f-K JO 

with 5*o = V%so- The latter condition means that the areas of the integration regions over s and s' in the two 
approaches are the same (see [^0] and [^6|,^lJ for more details). 

Using ([n]) and (|65| ) we can easily calculate the spectral density po(S, Q 2 ) 

9x S' 2 (25 + 3Q 2 ) , 

paiS, Q ) = — i 1) 67 

2(25 + Q 2 ) 3 

producing 

F™(Q 2 ) = -= s . (68) 

As one can see from the left part of Fig. || the difference between the two models in the region of interest (Q 2 > 
1 GeV 2 ) is very small. 



13 



0.5 



0.4 



. 3 



0.2 



. 1 



Q 2 F 7r {Q 2 ) [GeV 2 ] 

■ iiL t 




r 







10 



0.5 



0.4 



0.3 



0.2 



. 1 



Q 2 F„(Q 2 ) [GeV 2 ] 

,fl t 











10 



Q 2 [GeV 2 ] 



g 2 [GeV 2 



FIG. 5. Left part: Comparison of different LD models for the pion form factor. The solid line represents the tri angle LP 
model and the dashed line the square LD model. Right part: Comparison with available experimental data from fr3j73 ] of 
the pion form factor with (dashed line) and without (solid line) the Sudakov exponential in the 'triangle LD approach'. 



A. Sudakov effects due to the electromagnetic vertex in F 7 



TrLD (Q2) 



The crucial feature of the soft contribution is that it is accompanied by the Sudakov form factor. In other words, 
the double logarithms a s ln 2 (Q 2 ) do not cancel in this case. 

The one-loop radiative corrections to the spectral density p~o(S, t = Q 2 ) has been calculated by one of us (A.P.B.) 
J74| To analyze the Sudakov effects in the Feynman gauge, we need only the result for the gluon correction to the 
electromagnetic vertex (accompanied by the appropriate 1/2-insertions of self-energies into the quark lines, which 
gives an UV- finite result): 



vertex 



p(S,t) 



(a s C F /2TT)p (S,t) 



= 2 Li 



^)-»( i+ l>K 



i 



t 2 {6S + 5t) 
S 2 (25 + 3t) 



M2 + - 



1 2t {2S + 5t) 

2 + S (2S + 3t) 



(69) 



2 Li 2 



t 



2S + tJ ' t J'"" ') S 2 {2S + 3t) 

The leading asymptotics of this expression in the large t regime ispM 



t ) 8S 2 {2S + 3t) 



A EM-vertex^ 5; q 



po(S,t) t^oo 2tt 
So, we model the Sudakov corrections in the following way 



^ !l +2^ 



Olln- 



-Sudakov ( ^ g2) ^ - o( ^ Q 2 )exp 



oc s Cf 
2ir 



The modified spectral density is then used to model the soft term corrected by the Sudakov effects 



F 



TrLD-Sudakov 



(Q 2 ) 



l 



* 2 tt Jo 



So 



-Sudakov ( ^Q 2) d g 



(70) 



(71) 



(72) 



On the right part of Fig. | we show for comparison predictions for Q 2 Fj lLD (Q 2 ) and for Q 2 Fj rLD ~ Sudakov {Q 2 ) . 
One can see from this figure that the Sudakov effects in the electromagnetic vertex reduce (as expected) the soft 
contribution in the spacelike region by 6-20%. 



§ A similar correction was obtained in the light-cone QCD sum rule approach |J75l 
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B. Model dependence of the soft term in timelike region 



As we have seen in the previous subsection for spacelike values of the momentum transfer (Q 2 > 0) both LD 
models give rather close results for the pion form factor at Q 2 > 1 GeV 2 . But if we analytically continue these 
two models into the timelike region (q 2 = —Q 2 > 0), we obtain absolutely different results for both Re[F 7r (q 2 )] and 
Im^g 2 )]: 

Rc[F^(q 2 )} = 1 - 6(q 2 - As ) 1 " 6V ' 2 3/2 . (73) 

(i - 4s /r) 

Im[^D(? 2 )] = 0(4*0 - <? 2 ) 1 ~l 80/ t /2 , (74) 

(4s /Y - 1) 

Re[F™(g 2 )] = — 5 , (75) 

1 V2(l-q 2 /2S ) 2 

lm[Fj rLD (q 2 )} = . (76) 
We see that in the resonance region (q 2 < 4 GeV 2 ) the differ ences between these two models are rather large, and 




FIG. 6. Comparison of the analytic continuation to the timelike region q 2 — —Q 2 > for two different 
LD models of the pion form factor. The solid line shows the result for Jq 2 ^ (q 2 ) | , obtained from the 'square 
LD model', together with the real part of it (dotted line), whereas the prediction of the 'triangle LD model' 
for q 2 F 7r (q 2 ) is represented by the dashed line. The experimental data are taken from 

we can actually say nothing about the true behavior of F v (q 2 ) in this region. On the other hand, in the region 
q 2 > 6 GeV 2 the differences between the two models are less than the experimental uncertainties, and hence we 
can use them, at least as a first approximation, to model -F^g 2 ). Furthermore, in the case of the 'triangle LD' we 
have an explicit analytic expression for the Sudakov effects which we can now continue into the timelike region 



p^ e j^TrLD- Sudakov ^2 _ 
j m j^TrLD- Sudakov (g 2 )] = 



1 



1 



Sn 



& k ° v (S,q 2 ) dS; 



Sn 



-Sudakov 
"o,Im 



(5, q 2 ) dS , 



with 



-Sudakov / e „2 \ 
Po,Re \ h > 1 ) 



pX ak ° V (S, ^ = Ms ^ ^ sin i -a a tf)C F In 



p (S,q z )cos & s {q^)C F \n 



2S 



- 1 



Po(S, q 2 ) = p (S, -q 2 ) exp 



a s (g 2 )C F 

2ir 



q_ 

2S 
In 2 



1 

£_ 

2S 



(77) 
(78) 

(79) 
(80) 
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(Here d s (q 2 ) is the lowest-order model expression ( |l8| ) for a s in the timelike regime.) 

Using this model we obtain results, depicted on the lhs of Fig. |7[ After adding the analytically continued expres- 




FIG. 7. The lhs shows the analytically continued expression for the pion form factor to the timelike region q 2 = — Q 2 > 
with (dashed line) and without (solid line) the Sudakov exponential in the 'triangle LD approach'. The rhs shows the results 
for the total pion form factor in the timelike domain comprising the soft part within the 'triangle LD approach' and the 
hard one, calculated in [jrf^. Experimental data are taken from Jm . 

sion for the hard scattering (perturbative) part, including also transverse momentum effects (Sudakov+intrinsic 
effects) [ fT6[ , we arrive at the result, shown on the rhs of Fig. [?]. 



VIII. CONCLUSIONS 



In this paper, we investigated various aspects of the analytic continuation procedure from the spacelike to the 
timelike region of momentum transfers for several processes in QCD. We concentrated on studying several types 
of logarithmic contributions ln(Q 2 ) capable of producing ±Z7T in the timelike region. In the case of the ultraviolet 
logarithms, we reviewed the construction of the effective QCD coupling constant for the timelike region. The major 
result here is that the transition from a spacelike to the mirror timelike momentum only decreases the effective 
coupling constant. Studying the collinear logarithms, we established that in this case each eigenfunction 4> n {x) of 
the evolution equation acquires a phase factor e lSn . The phase vanishes for the asymptotic wave function, and there 
are no changes in this most realistic situation. But even in the case of the Chernyak-Zhitnitsky wave function, the 
interference effects are very small and, again, they decrease rather than increase the timelike contribution compared 
to the spacelike one. In the case of the pion electromagnetic form factor, we emphasized that the 7r 2 terms which 
may appear in the timelike region on the diagram by diagram level cancel in the total sum together with the double 
logarithms which generated them. Thus, we found no sources for the if-factor-type enhancements in the hard gluon 
exchange perturbative QCD contributions to the hadronic form factors. However, the situation completely changes 
if one considers the soft contribution. We investigated the simplest case of the pion electromagnetic form factor. 
To this end, we incorporated the local duality model suggested by the QCD sum rule studies performed earlier 
in the lowest (zero) order in a s . We included the a s correction which, as expected, contains the Sudakov double 
logarithms. In the timelike region the latter produce non-canceling positive it 2 terms which result in a if-factor- 
type enhancement. Our results for the soft contribution are in good agreement with existing experimental data on 
the pion electromagnetic form factor both in the spacelike and the timelike regions. We regard this agreement as 
another indication that soft contributions dominate the form factors at currently accessible momentum transfers. 
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